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SUMMARY 


A three-dimensional finite element solution of the acoustic vibration prob- 
lem in a solid propellant rocket motor is presented. The solution yields the 
natural circular frequencies of vibration and the corresponding acoustic pres- 
sure mode shapes, considering the coupled response of the propellant grain to 
the acoustic oscillations occurring in the motor cavity. The near incompressi- 
bility of the solid propellant is taken into account in the formulation. A 
relatively simple example problem is solved in order to illustrate the applica- 
bility of the analysis and the developed computer code. 


INTRODUCTION 


Present solid propellant rocket combustion instability state-of-the-art 
technology employs a linear analysis for predicting unsteady motions in the com- 
bustion chamber of the rocket based on the hypothesis that the influences of 
combustion and flow can be represented as perturbations of an acoustic problem 
in a closed chamber. In order to perform an accurate combustion instability 
analysis, one must have adequate knowledge of the natural circular frequency 
and corresponding pressure mode shape for the different modes of acoustic oscil- 
lation occurring in the motor cavity in the absence of combustion and flow. 

The cavities in functional rocket motors generally have complex shapes 
that may include a number of symmetrically located slots (or fins) cut into the 
propellant, the purpose of which is to increase the area of the burning surface 
(fig. 1). During the past several years, the need for experimental acoustic 
modeling of such irregular cavity geometries has been greatly reduced through 
the use of finite element modeling techniques (ref. 1). In the finite element 
method, a continuum is modeled as an assemblage of elements, connected at dis- 
crete locations (nodes) . Natural frequencies and mode shapes are then obtained 
from an eigensolution of the resulting equations of equilibrium of the modeled 
dynamic system. Applications of the finite element method to more generalized 
acoustic analysis problems are found in references 2 and 3. 

The presence of the solid propellant grain can significantly shift the 
acoustic system frequency from that of the gas phase alone, a portion of the 
acoustic energy being dissipated by the deformable solid material. This effect 
can be one of the more significant sinks for acoustic energy in both large and 
small rocket motors, the amount of damping depending on the grain geometry, pro- 
pellant mechanical properties, and the acoustic mode shape and frequency. It 
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is thus important to have the capability of finite element modeling of both the 
acoustic cavity and the propellant grain structure and, therefore, of a combi- 
nation and coupling of fluid and solid elements. Analyses of coupled plate- 
acoustic systems are found in references 4 and 5, and a three-dimensional finite 
element coupled acousto-mechanical analysis is formulated in reference 6. 

In the case at hand, the problem involves the coupling of the irrotational 
motions of an inviscid compressible fluid with the deformations of an incom- 
pressible (viscoelastic) solid material enclosing the fluid. Since the propel- 
lant grain is only accurately modeled as a nearly incompressible material, the 
well-known minimum potential energy finite element formulation, based upon 
Navier's equations of equilibrium and the Ritz procedure, is invalid. Formula- 
tions for cases of incompressible and nearly incompressible materials are found 
in references 7, 8, 9, 10, and 11. The solid finite element formulation utilized 
in this development is a linear displacement-linear mean pressure tetrahedron 
(ref. 10), similar to the Herrmann variational formulation (ref. 7) which employs 
a linear displacement function and a constant mean pressure function. 

The eigensolution technique employed in the developed program is that of 
condensation, which is discussed in reference 12. The technique utilizes the 
master and slave degree-of-freedom (d.o.f.) concept. In this application, the 
fluid pressure d.o.f. are designated master and the solid displacement and mean 
pressure d.o.f. are designated slave. 


SYMBOLS 


Values are given in both SI and U.S. Customary Units. The calculations 
were made in U.S. Customary Units. 
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exponential factors 
bulk modulus 

degree-of-freedom to strain transformation matrix 
constitutive coefficient matrix 

fluid element potential energy coefficient matrix 

fluid system potential energy coefficient matrix 

shear modulus 

mean pressure parameter 

identity matrix 

solid element stiffness matrix 

solid system stiffness matrix 

natural or volume tetrahedral coordinate 

solid element consistent mass matrix 

solid system consistent mass matrix 

pressure 

point (defined in fig. 2) 
surface area 

fluid element kinetic energy coefficient matrix 
fluid system kinetic energy coefficient matrix 
displacement components 
acoustic-solid coupling matrix 
element volume 

rectangular Cartesian coordinates 
direction cosine 
displacement vector 
strain 
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V Poisson’s ratio 

[$] transformation matrix defined by equation (14) 

p mass density 

0 ) natural circular frequency 


Subscripts: 


C 

J » 

p,q 

r,s 

S 


coupled 

element node identifiers 
mamber of degrees-of-freedom 
node identifiers (eq. (4)) 
solid 


Superscripts: 


e element 

T transpose (of a matrix) 


EQUATIONS OF EQUILIBRIUM 


Acoustic Cavity 

In applying the finite element method, the region of concern is divided 
into a number of subregions, or elements. Applications of the method in model- 
ing the behavior of a fluid continuum involve several pecularities not encoun- 
tered in the usual structural applications. In the discretization of a fluid 
continuum, the finite elements represent spatial rather than material subregions 
of the continuum; i.e., instead of representing finite elements of the fluid ma- 
terial, the elements represent subregions in the space through which the fluid 
moves (Eulerian description of motion) . Values of pressure at the nodal points 
of the element represent the pressures at the nodes rather than of the nodes 
(ref. 13). 

For a three-dimensional region, a volume element is associated with a num- 
ber of nodal points which define its shape. Within each fluid tetrahedral ele- 
ment the variation of the pressure p is prescribed by the values associated 
with the four element nodes through the expression: 




P(x,y,z) = ’"j j 


P5, 


( 1 ) 


where the natural or volume coordinates L^,L^,L^, and are defined mathemati- 
cally by 
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( 2 ) 


and physically by figure 1. 
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( 3 ) 


The matrix form of the homogeneous acoustic wave equations for the dis- 
cretized fluid contlnuiun Is given by 

([f] - (0^[T]){p} » {0} 

where F Is the coefficient matrix derived from a consideration of the potential 
energy of a fluid element, and X Is the coefficient matrix derived from a con- 
sideration of the kinetic energy of the fluid element. The natural circular 
frequency of the discretely modeled fluid continuum Is designated as b). 

The matrix F Is developed through a superposltlonlng of the Individual Sym- 
metric 4x4 element matrices, f, with restraint of the appropriate d.o.f. A typ- 
ical term of the f matrix Is (ref. 14) : 



, 9L 9L 

X X 





(4) 


where r and s take the values l,j,Jt,m cyclically. The fluid mass density and the 
volume of the tetrahedron are designated p and V, respectively. 


The matrix T Is likewise developed through the superposltlonlng of the In- 
dividual symmetric 4x4 element matrices, t, with restraint of the same appro- 
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where B is the bulk modulus of the fluid, 
evaluating the terms in the t matrix is: 


A useful integration formula for 


/ ,a,b^c^d ^ ^ j a!b!c!d! 

Vj Vm ■ (a+b+c-M+3) ! 


( 6 ) 


Solid Propellant 


Within each solid tetrahedral element, the variation of displacement is 
prescribed by the displacement of each element node through the expression: 


u (x , y , z) 
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( 7 ) 
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where I is the 3x3 identity matrix and where 


L'^iJ * L“i ''i "ij’ 


( 8 ) 


The u,v,w terms are the displacement components in the x,y,z directions, respec- 
tively. 

The matrix form of the homogeneous equations of motion for the discretized 
solid continuum is given by 


( 




(9) 


where K is the stiffness matrix and M is the mass matrix- The term 0 )g is the 
natural circular frequency of vibration of the solid system. The linear varia- 
tion of the mean pressure parameter 


H(x,y,z) = 


-z — :p— (e + e + e ) 
l-2v XX yy zz 


( 10 ) 


is prescribed in the same manner as were the pressure (in the acoustic cavity) 

‘h. 


and the displacement function, i.e*, 

H(x,y.z) = = L^i L. LJ 
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( 11 ) 


The terms e , e , and e are the normal strain components and V is Poisson’s 
xx’ yy* zz 

ratio. 

The matrix K is developed through a superpositioning of the individual 
symmetric 16x16 element stiffness matrices, k, with restraint of the appropri- 
ate d.o.f. The form of the matrix k is given by (ref. 10): 


[k] = 


D^j 0 


0 


[$]'*'[e][$] dxdydz 


D [ 0 

n I T 


( 12 ) 


where I is the 4x4 identity matrix and where 
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The matrix M Is developed through the superpositioning of the individual 
symmetric 12 x 12 element mass matrices, m, with restraint of the same appropri- 
ate d.o.f. The consistent mass matrix, derived in references 15 and 16, is 
used in the formulation and is defined by: 

[m] * Pg [l]^[l] dxdydz 

where Pg is the mass density of the solid propellant grain and where 

[l] = [L^I Ljl L^I L^l] 

Acoustic-Solid Coupling 

The formulation of the coupled acoustic-solid equations through equilib- 
rium considerations is foiond in references 5 and 6 . Reference 5 deals with the 
transient response of a coupled system while, in reference 6 the eigenvalue 
problem is formulated. The matrix form of the homogeneous equations of equi- 
librium for the discretized coupled system is given by 
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(18) 


where is the natural circular frequency of the coupled system. 

The matrix U of equation (18) has dimensions pxq, where p is the number of 
pressure d.o.f. and q is the number of displacement and mean pressure parameter 
d.o.f. The components of the U matrix are projections of the cavity-propellant 
interface area associated with each surface nodal point, the surface area dis- 
tribution being physically defined in figure 3. The portion of the U matrix 
associated with the i,j,il triangular surface area of a boundary tetrahedron is 
given by 
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where ••• 9 direction cosines and S^, S, and are the con- 
tributing (to a node) surface areas. ^ ^ 


SOLUTION 

Solution of the system mathematically modeled by equation (18) is facili- 
tated by the use of condensation (ref. 12). In applying the technique to this 
problem, the pressure d.o.f. are designated master and the displacement and 
mean pressure parameter d.o.f. are designated slave . 
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Example Application 


The developed program is applied to the solution of the simple acoustic- 
solid system shown in figure 4. The plane defined by nodes 5678 is the 
acoustic-solid interface. The pressure is set equal to zero at nodes 1,2,3, 
and 4. Thus, for the corresponding solution, p = 4 and q = 32. ^ 

The properties of the cavity gas used are: a mass density p of 1.225 kg/m 

(1.146 X 10~^ Ib-sec^/in^) and a bulk modulus B of 1.42 x 10^ Pa (20.59 Ib/in^). 
The properties of the solid propellant grain used are: a shear modulus G of 

7 2 

4.41 X 10 Pa (6400 Ib/in ), a Poisson^s ratio V of 0.5 and a mass density Pg 

of 1.770 X 10^ kg/m^ (1.656 x lO"^ lb-sec^ / in^) . 

The solution for the coupled system yielded frequencies of 488 Hz, 1008 Hz, 
1193 Hz, and 1804 Hz, compared to respectively corresponding values of 3670 Hz, 
6934 Hz, 10,358 Hz, and 12,368 Hz obtained for the fluid system alone. 


CONCLUDING REMARKS 

It has been demonstrated through the use of a simple example that the 
developed computer code can be used to predict the effect of solid propellant 
grain structure on the natural frequency associated with the acoustic vibra- 
tions occurring in the rocket motor cavity. It can be concluded that the effect 
of certain suppression devices such as resonance rods can be predicted in the 
same manner. 
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Figure 1.- Solid propellant configuration. 
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Figure 3.- Surface areas associated with 
interfacial nodal points. 



Figure 4.- Acoustic-solid system. 
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